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Quantum ladder models, consisting of coupled chains, form intriguing systems bridging one and two dimen- 
sions and have been well studied in the context of quantum magnets and fermionic systems. Here we consider 
ladder systems made of more exotic quantum mechanical degrees of freedom, so-called non-Abelian anyons, 
which can be thought of as certain quantum deformations of ordinary SU(2) spins. Such non-Abelian anyons 
occur as quasiparticle excitations in topological quantum fluids, including px + ipy superconductors, certain 
fractional quantum Hall states, and rotating Bose-Einstein condensates. Here we use a combination of exact 
diagonalization and conformal field theory to determine the phase diagrams of ladders with up to four chains. 
We discuss how phenomena familiar from ordinary SU(2) spin ladders are generalized in their anyonic coun- 
terparts, such as gapless and gapped phases, odd/even effects with the ladder width, and elementary 'magnon' 
excitations. Other features are entirely due to the topological nature of the anyonic degrees of freedom. In 
general, two-dimensional systems of interacting localized non-Abelian anyons are anyonic generalizations of 
two-dimensional quantum magnets. 

PACS numbers: 74.20.Mn, 67.80.kb, 75.10.Jm, 74.75.Dw, 74.20.Rp 



I. INTRODUCTION 

Quantum antiferromagnets and, more generally, electronic 
systems are notoriously known to behave fundamentally dif- 
ferently in one and two spatial dimensions. In one spatial 
dimension quantum fluctuations are enhanced and often give 
rise to critical properties such as algebraic spin (or charge) 
correlations^. In contrast, in two spatial dimensions one fre- 
quently finds quantum ground states with long-range order, 
which often originate from the spontaneous breaking of a con- 
tinuous symmetry. The archetypal example of the latter are 
Heisenberg antiferromagnets on bipartite lattices^, where the 
Neel ground state arises from the spontaneous breaking of 
the SU(2) spin synmietry and the resulting zero-energy Gold- 
strone modes are spin- wave excitations, also called magnons. 
Quantum ladder systems, consisting of a finite number W 
of coupled one-dimensional (ID) systems, form a bridge be- 
tween between these two limits, and the evolution of quan- 
tum ground states in the dimensional crossover of increasing 
ladder width has been well studied in the context of itiner- 
ant bosonic and fermionic systems^ as well as quantum spin 
ladderP. A variety of remarkable crossover effects have been 
observed, such as the celebrated even/odd effect in quantum 
spin- 1/2 Heisenberg ladders, where gapless ground states are 
found for all odd width W, while ladder systems with an even 
number of legs exhibit a spin gap"^. Surely, this effect has 
also been experimentally observed for actual materials realiz- 
ing almost perfect two and three leg S=l/2 antiferromagnetic 
(AFM) ladders when performing careful magnetic suscepti- 
bility measurements^ . 

In this manuscript, we consider systems of more ex- 
otic quantum mechanical degrees of freedom, so-called non- 
Abelian anyons, which have attracted considerable interest 
in the description of non-Abelian vortices in unconventional 
Px + Wy superconductors^, quasiholes in certain fractional 



quantum Hall stateJ^, or vortices in rotating Bose-Einstein 
condensates^^, and in theoretical proposals for inherently fault 
tolerant quantum computing schemes Since we are inter- 
ested in their collective quantum ground states in two spatial 
dimensions, we follow a route similar to the above-mentioned 
studies of SU(2) quantum antiferromagnets and study ladder 
systems of interacting anyonic degrees of freedom. 

Formally, non-Abelian anyons can be described by so- 
called su(2)/g Chern- Simons theories, which correspond to 
certain quantum deformations^^ of SU(2). In these theories, 
the non-Abelian degrees of freedoms are captured by 'gener- 
alized angular momenta' j, which for a given su(2)k theory, 
are constrained to the first k -\- 1 representations of SU(2) 

. _ 1 3 k 
~ ' 2' ' 2'*"' 2 ' 

Similar to the coupling of ordinary angular momenta, two 
non-Abelian degrees of freedom can be 'fused' into multiple 
states with total angular momenta (or spins) 

ji^j2 = |ji-j2|©|ji-j2| + l©. . .©min(ji+j2,/^-ji-j2) , 

where again the 'cutoff k of the deformation enters. The oc- 
currence of multiple fusion channels on the right-hand- side of 
the above equation is what intrinsically gives rise to a macro- 
scopic degeneracy of states for a set of multiple non-Abelian 
anyons - the hallmark of non-Abelian statistics. 

In this manuscript we consider the fundamental case of 
non-Abelian anyons with generalized angular momentum j = 
1/2, which obey the fusion rule l/2(g)l/2 = 0©1 remi- 
niscent of two ordinary spin-l/2's coupling into a singlet or 
triplet. Like the Heisenberg Hamiltonian for ordinary spins, 
interactions between the anyons energetically split the two fu- 
sion outcomes, which in the case of ordinary SU(2) spins is 
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captured by the familiar Heisenberg Hamiltonian 

H = ^^^JijSi-Sj (1) 

= Ji^n^^^. + const., (3) 

which can be viewed as a sum of pairwise projectors 11^ ^ onto 
the singlet state. Similarly, we can define an 'anyonic Heisen- 
berg Hamiltonian' that for a pair of non-Abelian anyons with 
generalized angular momentum j = 1/2 projects onto the 
j = (singlet) fusion channel, thus taking the same form 
as Eq. ([3]) of the Hamiltonian above. In analogy to ordinary 
SU(2) spins we refer to positive couplings (projecting onto the 
generalized j = state) as 'antiferromagnetic', while nega- 
tive couplings (projecting onto the generalized j = I state) 
are 'ferromagnetic'. 

It has recently been shown that similar to their SU(2) coun- 
terparts chains of interacting non-Abelian anyons can exhibit 
a variety of collective ground states including stable gapless 
phase^SEH and exotic infinite-randomness fixed points 
In this manuscript, we aim at understanding two-dimensional 
ground states of interacting non-Abelian anyons and - follow- 
ing a similar route as in the case of the above-mentioned stud- 
ies of SU(2) quantum antiferromagnets - we consider systems 
of coupled chains forming VK-leg ladders. Employing exten- 
sive numerical simulations combined with a conformal field 
theory analysis, we investigate phase diagrams of I^-leg lad- 
der with up to = 4 legs, which allows us to derive some 
conclusions also for the 2D limit of W oo. We mostly 
focus on the case of su(2)k with /c = 3 as a representative 
example and in particular all numerical simulations are per- 
formed for this case. We will return to the more general case 
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of arbitrary level /c > 3 in Sec. VI 



While in this manuscript we detail the physics of interacting 
non-Abelian anyons mostly in terms of (deformed) quantum 
spins - a notion more familiar to the field of low-dimensional 
quantum magnetism, we have put forward another perspective 
on the physics of interacting anyons in the context of certain 
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FIG. 2: (color online) Representation of the parameter space of the 
model on a circle. The couplings Jrung and Jieg can be either pos- 
itive (AFM) or negative (FM). The different phases in each of the 
quadrants labelled by " Jrung "-"-^leg" are given in Table|l| The (grey) 
squares at 7r/4 correspond to isotropic couplings | Jrung /Jieg| = 1. 
The (orange) stars mark the parameters in the vicinity of the strong 
coupling limit | Jrung| ^ |Jieg| used in this work, namely Q — 
{37r/7, 47r/7, IOtt/T, IItt/T} labeled from 1 to 4, respectively. 



fractional quantum Hall states in a recent articlPI. There we 
have made a connection between the collective states of (any- 
onic) excitations in non-Abelian quantum Hall liquids and the 
physics of moving on a non-Abelian quantum Hall plateau. 

The remainder of this manuscript is structured as follows: 
We will start with a detailed derivation of the microscopic 
models analyzed in this manuscript in Sec.[ll| This is followed 
by a discussion of the phase diagrams of the various ladder 
models starting from the strong rung-coupling limit in Sec.[lll| 



and continuing with the weak rung-coupling limit in Sec. IV 
We will then turn to the peculiar role of boundary conditions 
and the occurrence of gapless modes at open boundaries for 
these anyonic ladder models in Sec. [V] We round off the 
manuscript by a discussion of the two-dimensional limit of 
these ladder models in Sec . |VII| and generalization to su(2)fc 
theories with /c > 3 in Sec.lVIl 



II. THE MICROSCOPIC LADDER MODEL 



Jleg 




FIG. 1 : Three coupled chains of interacting anyons (indicated by the 
filled circles). The interaction along ( Jieg) and perpendicular ( Jmng) 
to the chains are indicated by the ellipses. 



In this section we will give a definition of the microscopic 
W^-leg ladder models for so-called su(2)3 Fibonacci anyons. 
We will keep our discussion short but self-contained, as an 
extended derivation of general microscopic Hamiltonians has 
been given in Ref. 21 We will emphasize in the following 
those aspects that are not covered in Ref. 21 For a given 
VK-leg ladder we denote the strength of the interactions as 
Jleg and Jrung for the coupling along and perpendicular to 
the chains, respectively, as illustrated in Fig.[T] Parametrizing 
these couplings as Jieg = cos d and Jrung = sin d we will map 
out the parameter space on a unit circle as shown in Fig. [2] 

Our numerical analysis of these ladder systems is based on 
exact diagonalization using the Lanczos algorithm which pro- 
vides us with the low-energy spectra of finite systems with 
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extent W x L, where L is the length of the ladder in the chain 
direction and W is the width of the ladder in the rung direc- 
tion. In our exact diagonaliztion studies, we have been able to 
analyze systems of size 2 x L (L = 8, 12, 15, 16, 18, 20, 21), 
3 X L (L = 6, 8, 9, 10, 12) and 4 x L (L = 4, 6, 8). 
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A. The basis states 

To describe the basis states of a set of N localized (inter- 
acting) su(2)3 anyons we consider a fusion path, as shown in 
Fig. [3ja). The basis of the many-anyon Hilbert space cor- 
responds to all admissible labelings ^2, . . .) of the links 
in this fusion path with labels Xi corresponding to general- 
ized angular momenta of su(2)3. These labelings must satisfy 
the constraints of the fusion rules at each vertex of this fusion 
path. For su(2)3 these fusion rules are 
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where a g) /3 = /3 g) a. These fusion rules ^ reveal an au- 
tomorphism a^d = 3/2 — a, allowing an identification of 
0^3/2 and 1/2^1 for su(2)3. Using the notation for the 
Fibonacci theory, we write the identity 1 for the former and 
the label r for the latter, thus leading to the fusion rules 



lg)l = 1 

1 g) r = r 
r g) r = 1 



(5) 
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For the labelings of the fusion path these rules then imply that 
Xi = 1 has to be followed by x^+i = r but Xi = r can be fol- 
lowed by either x^+i = 1 or x^+i = r. This constraint gives 
an overall Hilbert space size to Fat+i -\-Fn-i ~ (j)^ (for large 
N) where F is the Fibonacci sequence and = (1 + a/5)/2, 
the golden mean. Note that in comparison to ordinary SU(2) 
spin- 1/2 systems this Hilbert space has a reduced size. 

Our specific choice of a fusion path for the VK-leg ladder 
system is shown in Fig. |4] Using periodic boundary condi- 
tions along the leg direction enables us to conveniently use 
the translation symmetry of the system along the legs: The 
Hamiltonian matrix can then be block-diagonalized into L 
blocks labeled by the total momentum = 27r ^ of the eigen- 
states. Hence, the Hilbert space (in each symmetry sector) 
grows approximately as (j)^ / L which is one of the limiting 
factors of our simulations. To provide some examples, the 
Hilbert spaces of the i^T = sector for 2 x 21, 3 x 12 and 
4x8 ladders are found to be of sizes 28 527448, 2 782 659 
and 609 147, respectively. 



FIG. 3: (color online) a) the fusion path connecting N Fibonacci 
(su(2)3) anyons. Basis states correspond to all admissible labelings 
of the edges Xi with anyon charges and 1/2 that satisfy the fu- 
sion rules, b) a nearest neighbor coupling in a chain of anyons can 
be calculated by using the F matrix to transform to a new basis in 
which the fusion product of the two anyons is one of the variables, c) 
longer-range interactions first need to be mapped to nearest neighbor 
interaction by braiding anyons. d) a Dehn twist, giving an additional 
phase factor is needed if the interaction winds around the torus. For 
su(2)3 the phase is V^cc = for x = and '0cc = — 47r/5 for x = 1/2. 



B. The rung interactions 



With our choice of fusion path the rung coupling Jrung on 
ladders with open boundary conditions on the rungs always 
connects neighboring anyons along the fusion path. To cal- 
culate the interaction between two neighboring anyons as in 
Fig.[3];b) we need to calculate their total spin by performing a 
basis transformation using the F-matrix, and then assign en- 
ergy — Jrung to the identity fusion channel and energy to the 
r fusion channel. This basis transformation is illustrated in 
Figure [3] 

Denoting the local basis states on the three edges around the 
interaction as x^, x^+i) G {|l,r, 1), |1, r,r), |r,r, 1), 

|r, l,r), |r, r, r)} and the states after the F-transformation 

as \xi-i,Zi,XiJ^i) e {|1,1,1), |1, r,r), |r,r,l), |r,l,r), 
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FIG. 4: (color online) Sketch of a typical ladder system of extent 
4 X L (L = 6), where the dots indicate the location of r-anyons. The 
fusion path C (dotted line) with labelings Xi (or yi) is used define an 
(arbitrary) ordering of the sites, which is used in the definition of the 
basis states. The exchange couplings Jieg and Jrung are indicated. 



= r the local Hamiltonian is Hi = —FiPiFi where Pi is 
the projector onto the state with = 1. In the basis defined 
above we get 
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The rung Hamiltonian is obtained by multiplying this ma- 
trix by Jrung and, for the term shown in Fig. [4] acts on the 
local states 1^1,^2,^3). 



C. The leg interactions 



|r, r, r)} we can write the F-matrix as 
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Assigning an energy —1 to the identitiy 



The leg couplings shown on Fig. |4j on the other hand, are 
longer range interactions and requires 'braiding' of anyons as 
illustrated in Fig. |3jc). Let us first consider a next-nearest 
neighbor interaction along a chain. To transform this into 
(6) a nearest neighbor interaction we need to change the basis 
again, this time by braiding the two left anyons in a clock- 
wise manner with a so-called braid matrix Bi acting on the 
states Xi+i). Using the same basis as before this 

1 and to braid matrix can be written as: 
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(8) 



and the next nearest neighbor coupling then becomes 

B\Hij^iBi. 

Similarly, for the leg coupling, illustrated for a four-leg lad- 
der in Fig. [4] we need three braids and act on the whole se- 
quence X2, X3, X4, X5, xe) involving 6 bonds (in general, 
involving 2W — 2 bonds for W chains) along the fusion path 
C according to the linear transformation: 

^leg^ = ^leg B2^B^^B^®n^^Bl® bI bI . (9) 

Similar formulas can easily be derived for any bond and any 
width W. 

It should be noticed that acting on any given initial state 
\xiX2 - ■ ■ xlw)^ each of i\\tW • L leg couplings can poten- 
tially generate up to 2^^~^ resulting linear independent states 
since each operator in ([9]) can generate up to two such states. 
This exponentially growing number of resulting states should 
be contrasted to the single state generated by a spin flip oper- 
ation in the case of ordinary SU(2) spins. As a consequence. 



this leads to denser and denser matrices for increasing W in 
the anyonic ladder models, which limits the numerically ac- 
cessible system sizes for larger width W . 



D. Periodic boundary conditions along the rungs 

Closing these open boundaries along the rung direction is 
done by adding additional couplings between the first and last 
legs as shown in Fig. |5] To calculate the Hamiltonian ma- 
trix for these couplings one first has to again braid the two 
involved anyons until they are nearest neighbors along the fu- 
sion path. The subtlety with this term is that after the braid- 
ings we do not end up with the usual nearest neighbor term 
of Fig. |3jb) but with a coupling that twists once around the 
fusion path as illustrated in Fig. [3jd). Untwisting this wind- 
ing by a 27r rotation of the right anyon and all following ones 
by 27r around the fusion path gives rise to a Dehn twist phase 



5 



1.6 



If 



rung 



;F4 



c 



c 



-•-r 9-^ •- 



FIG. 5: (color online) A ladder system as in Fig.|4] but with peri- 
odic boundary conditions, where an an additional rung coupling term 
connects the two out legs (vertical red segments). 



factor ex.p{i'^xj^, which is 1 for x = 1 but exp(— 4z7r/5) for 
X = r. 

The Hamitlonian for this rung term acts on the local sites 
^1,^2,^3,^4,^5 of Fig. [5] and reads 

^ru"n| = ^mng ^2 ^ ^3 ^ ^4 ^ T] ^ BI ^ BI , (10) 



where Ti is the (diagonal) 2x2 twist matrix, 
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e-"^*^/^ 



(11) 



in the local {|1), |r)} basis (for the variable Xi). 

Care must be taken in choosing a consistent convention for 
the phase of (counter)-clock wise braids and Dehn twists. An 
inconsistent choice can easily be detected as it will cause a 
broken translation symmetry along the rungs that can be seen 
in, e.g. the local bond energies. 
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FIG. 6: (color online) Dispersion of a generalized 'magnon' exci- 
tation along the ladder direction for two-leg and four-leg ladders 
and various couplings. The magnon excitation is created by flip- 
ping a local label Xi from r to 1, which can then propagate down 
the ladder. Data is shown for 2 x 20 (and 2 x 21) ladders (open 
symbols) as well as 4 x 8 ladder (closed symbols) and coupling 



parameters 6 — n/A (J\eg = Jn 



a/2) and 6* = 37r/4 



(Jieg = — Jrung = — y 1/2). Thc soUd Hncs are a guide-to-the- 
eye, obtained from Fourier series fit to the data. The arrows indicate 
extrapolations to the thermodynamic limit as shown in Fig.jS] 
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TABLE I: The various phases of the Fibonacci ladders as a function 
of the number of legs W. Each line corresponds to one of the four 
quadrants of the parameter space shown in Fig|2] In the first column, 
the first (second) label refers to the rung (leg) coupling Jrung (-^leg). 
Gapped phases are labeled by A. For the gapless phases the value of 
the central charge of the low-energy CFT is indicated. 



III. STRONG COUPLING LIMIT AND PHASE DIAGRAMS 

For SU(2) quantum spin ladders a single phase extends 
from the weak to strong rung coupling limit for any of the 
four possible signs of the rung and leg couplingPES. The 
generic phase diagram thus has at most four different phases. 
For su(2)3 anyonic ladder we observe the same behavior and 
we will start to discuss the various phases starting from the 
strong rung coupling limit | Jrung | ^ | J\eg \ • The results dis- 
cussed below are summarized in Table |l| The total spin of an 
isolated rung, which depends on the sign of the rung coupling 
and on the rung length W, completely determines the nature 
of the phase at finite Jieg and whether it is gapped or critical. 
For antiferromagnetic Jrung (first two lines of Table [I]) we find 
similar even/odd effects as in the SU(2) case. Even widths 
are gapped while odd widths are critical and characterized by 
the same CFT as the single chain. For ferromagnetic Jrung 
and W = 3p (p an integer) the rungs form singlets, (labeled 



with the identity 1) and hence, the ladders are gapped. Oth- 
erwise, the rungs behave as "triplet" (r) states and the low- 
energy physics is that of an (effective) critical chain as shown 
in the two last lines of Table Jl 



A. Antiferromagnetic rung coupling 

Let us first consider AFM rung coupling, where for an iso- 
lated rung the ground state has total angular momentum j = 
(state with label r) for even width W and j = ^ (state with 
label 1) for odd width. 

For even W, the ground state at Jieg = is a prod- 
uct of local r states on the individual rungs. The elemen- 
tary excitation is a local singlet (1) excitation with a gap 
Aq{W) ~ 1/W. For (weak) leg coupling the elementary 
'magnon' excitation can hop to one of its two neighboring 
rungs in first-order in Jieg, giving rise to a dispersion of 
width ex I Jieg I . Typical such dispersions (but for intermedi- 



6 



a)^=0 h)K=n 




0.04 0.08 0.12 0.16 0.04 0.08 0.12 0.16 
inverse system size 1/L inverse system size 1/L 



FIG. 7: (color online) Finite-size extrapolations of the energy gaps 
of 3-leg ladder systems for various couplings 0, with ferromagnetic 
rung and antiferromagnetic leg couplings in the sl) K = and 
h) K = TT momentum sectors. Extrapolations to the thermody- 
namic limit are obtained by fitting the numerical data to the form 
A(L) A((X)) + ^ exp (— L/^), where ^ is a correlation length. In 
the = sector, the first excitation energy (open symbols) extrap- 
olates to zero (indicating that the ground state is two-fold degenerate 
in the thermodynamic limit), while the extrapolation of the second 
excitation energy (filled symbols) indicated a finite gap. 
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FIG. 8: (color online) Extrapolated values A(L ^ oo) of the gaps 
at the two crystal momenta K = and K = Kq corresponding to 
the zero-energy modes of the single chains. A(oo) is plotted as a 
function of coupling parameter 0, i.e. G [0, tt] for 2- and 4-leg 
ladders with AFM rung coupling (a) and G [— vr, 0] for a 3-leg 
ladder with FM rung coupling (b). For AFM (FM) leg coupling, 
Ko = n (Ko = 27V / 3) and clusters up to 2 x 20 (2 x 21), 3 x 12 
(3 X 12) and 4x8 have been used. A(oo) is deduced by fitting the 
data as A(L) ~ A(oo) + ^ exp (— L/^), where ^ is a correlation 
length. 



ate couplings) are shown in Fig. [6] The gap decreases linearly 
as Ao(I^) — Q^l Jiegl + O(Jigg). However, this perturbative 
strong-coupling result for the gap is restricted to a shrinking 
region ~ 1/W around Jieg = as 1^ gets larger. For ordi- 
nary SU(2) ladders it has been argued that the gap vanishes 
as exp (—cW) for large enough W and any chain coupling 
Jieg We will return to the question whether a gap can sur- 
vive for anyonic systems in the limit 1^ ^ oo below. 

Away from the above-discussed limit, the gaps of the 2- 
and 4-leg ladders can be obtained for intermediate couplings 
by using Lanczos exact diagonalisations of clusters of differ- 
ent lengths. Finite size scalings (similar to the one shown in 
Fig. |7] for a 3-leg ladder to be discussed later) enable to ac- 
curately estimate, in the thermodynamic limit, the gaps at the 
minima of the dispersion (see e.g. Fig. |6]) of the excitation 
spectrum. Results of the extrapolated gaps are summarized 
in Fig. [SJa). Note that the minima of the dispersion occurs 
at different momenta depending on the sign of Jieg, and tt 
for antiferromagnetic Jieg > 0, and 27r/3 for ferromagnetic 
Jieg < 0. In the latter case, for sufficiently large leg coupling, 
the minima at 27r/3 can disappear as shown in Fig.|6] 

For odd width W, since the GS of a single AFM rung car- 
ries angular momentum j = 1/2 (r), the low-energy effective 
model for weakly coupled rungs is that of a single r-anyon 
chain. Indeed, as shown in Fig. |9] for a 3-leg ladder we find 
that the low-energy spectrum is gapless and can be described 
by a Conformal Field Theory (CFT) identical to the one of a 
single chain^^. In particular, we find that the lowest energies 



(per rung) scale as 

c 1 

en{L) eoo + ^^(~ ^2 ^ ^^^^^^ ' ^^^^ 

where L is the length of the ladder, uis 3. zero-mode velocity, 
and 2hn and c are the conformal weights (or scaling dimen- 
sions) and central charge of the CFT. Depending on the sign of 
the leg coupling these gapless theories are those of the tricrit- 
ical Ising model (c = 7/10) or 3-state Potts model (c = 4/5) 
for AFM and FM couplings respectively. 

The anyonic ladders with AFM rung coupling thus behave 
similarly to their SU(2) analogs with even/odd widths giving 
rise to gapped/gapless physics as summarized in the first two 
lines of Table [D 



B. Ferromagnetic rung coupling 

Next, we move to the case of a ferromagnetic rung cou- 
pling, where we find major differences between the anyonic 
ladders and their SU(2) counterparts. For ordinary SU(2) lad- 
ders of even width W, the strong coupling rungs form a total 
integer spin and the effective low-energy model (for weak leg 
coupling) is a Haldane Heisenberg chain^"^ which is gapped 
for AFM leg coupling. For odd width W the ladders remain 
gapless for either sign of the leg coupling, since each rung 
forms a r state. 
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a) 2-leg ladder, 4-leg ladder b) 3 -leg ladder 
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FIG. 9: (color online) Finite size scaling of 3 x L (3 -leg) ladders of 
sizes up to L = 12, with AFM rung coupling and with both AFM 
(top) and FM (bottom) leg couplings. The lowest eigen-energies per 
rung (symbols) are plotted vs The right and left panels corre- 

spond to two different crystal momenta, the ones characterizing the 
zero-energy modes of a single chain. The values (indicated on the 
plot) correspond to the strong rung couplings | Jrung/Jieg| ^ 4.38 
labelled as '1' and '2' in Fig.|2] CFT scalings are shown: (i) a linear 
fit of the GS energy (labeled by '0') accurately provides the overall 
energy scale (i.e. the velocity); (ii) all other straight lines are ex- 
pected CFT scalings using the conformal weights 2hn and the central 
charge c indicated on the plots (see Eq. |T2] ) in the text). 
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FIG. 10: (color online) Finite size scaling of 2 x L (2-leg) ladders 
with FM rung coupling, with both FM (top) and AFM (bottom) leg 
couplings and with up to 2 x 21 and 2 x 20 sites, respectively. The 
values (indicated on the plot) correspond to the strong rung couplings 
|-^rung/-^ieg| ^ 4.38 labelled as '3' and '4' inFig.[2] Same notations 
and same analysis as in Fig. [9] 





FIG. 1 1 : (color online) Phase diagrams of the 2-leg and 4-leg ladders 
(a) and of the 3-leg ladders (b) versus the couplings Jieg = cos and 



J ri 



sin 0. The central charge of the gapless phases is indicated. 



Also in contrast to ordinary SU(2) spins, for anyonic lad- 
ders, we find different phases and periodicities of 3 in for 
/c = 3. Ladders with = 3p (p an integer) and ferromag- 
netic Jrung are gapped, since each rung forms a singlet (1) 
state similar to the even width ladders in the AFM case. As 
an example, we show in Fig.jSjb) the gap of a 3-leg ladder ob- 
tained from finite size scalings, examples of which are shown 
in Fig. [7] Alternatively, the low-energy effective model of lad- 
ders with widths that are not multiples of 3 is again that of a 
single r-chain and thus gapless as illustrated in Fig.[TO|for a 

2- leg ladder. One might naively expect that the 2-leg ladder is 
again a gapped Haldane chain, since two FM coupled j = 1/2 
momenta form a total j = I momentum. However, as noted in 
the introduction, in su(2)k theories with odd level k one can 
identify momentum j with momentum k/2 — j by fusing it 
with the Abelian momentum-/c/2 particle. For su(2)3 this im- 
plies that momentum j = I behaves like momentum j = 1/2 
(r). We find that this gapless phase extends all the way up to 
weak rung coupling. 

We summarize our results in the phase diagrams of 
Fig. [TTJa) for 2-leg and 4-leg ladders, and of Fig. [TTJb) for 

3- leg ladders. 



IV. DECOUPLED CHAINS 

We now turn to a discussion of the limit where the rung 
coupling between the individual legs of the ladder vanishes. 
In contrast to the case of conventional SU(2) spin ladders, we 
find that the anyonic ladder system does not decompose into 
independent chains in this limit of vanishing rung coupling, 
i.e. Jrung = 0. In particular, we find that the energy spec- 
trum in this limit is not given by the free tensor product of 
the energy spectra of individual chains, but rather turns out 
to be a certain subset thereof. In the following, we describe 
a set of 'topological gluing conditions' that constrain the en- 
ergy spectrum to this subset of the free tensor product. We 
closely follow the analytical arguments, which we developed 
in Refs. ll5l20l in a so-called 'liquids picture', where we iden- 
tify the collective gapless modes of the quasi one-dimensional 
anyon chains (or ladders) with edge states at the spatial inter- 
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2 2 

1/L 1/L 

FIG. 12: (color online) Finite size scaling of the lowest eigen- 
energies (normalized per rung) of a 2-leg ladder with Jrung = and 
Jieg = ±1 (decoupled chains) vs . The right and left panels 
correspond to two different crystal momenta, the ones characteriz- 
ing the zero-energy modes of a single chain. Same procedure and 
notations as Fig.[9]and Fig.[TO| 



face between two distinct topological quantum liquids - for 
an illustration see e.g. Fig. 2 of Ref. [151 This 'liquids pic- 
ture' provides a set of analytical rules which allow to obtain 
the spectrum of these decoupled any on chains, which in the 
remainder of this section we compare with numerical results 
for 2-leg and 3 -leg ladder systems. We find perfect agreement 
of the two approaches. 

Let us briefly describe the analytical spectrum of the de- 
coupled chains (Jrung = limit), based on the results ob- 
tained in Ref. T3l(see Fig. 2 of that reference). At each in- 
terface between two topological (or Hall) liquids there is an 
edge (see Fig. |2jb) of Ref. 15 ). The key tool developed in 
Ref. [15] was that each chain can be viewed as 'filled' with 
a new topological (or Hall) liquid so that the right- and left- 
moving gapless degrees of freedom of each chain arise from 
the juxtaposition of two topological liquids. The field theory 
describing each of these edges arises from the familiar GKO 
CO set construction^^ of conformal field theory. Consider for 
example two decoupled chains. Thus, there are five liquids, 
and four edges. For, say, "antiferromagnetic (AF)" interac- 
tions between the anyons, the spectrum of these four edges in 
the 'topological sector' of 'topological charge' ji takes on the 
form (following the rules developed in Ref. tl5j and using the 
notation of the same article) 

{i^Lt ^R)i ^lyI (Vfi)jj 

Here ji denotes the 'topological charge' which is "ejected" 
from the four-edge system to infinity through the surround- 
ing ("parent") topological liquid [compare again Fig. |2jb) of 
Ref. J5J . The left (holomorphic) and right (antiholomorphic) 
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TABLE II: Scaling dimensions of for two (top) and three (bot- 
tom) decoupled chains with antiferromagnetic (0 = 0) and ferro- 
magnetic (0 = tt) leg couplings, for the two crystal momenta cor- 
responding to the respective zero-energy modes of a single chain. 
Listed here are the numerically observed conformal operators in the 
subset of the free tensor product of energy states for W individ- 
ual chains, corresponding to conformal field theories with central 
charge c = W x (c = VKx|)for antiferromagnetic (fer- 
romagnetic) coupling. For the Ising theory (c = 7/10), we use 
the common identification of operators with conformal weights, i.e. 
/ ^ 0,e ^ 1/5, ^ 6/5, e'' ^ 3, cr ^ 3/40, cr' ^ 7/8. 
For the 3-state Potts model (c = 4/5) this identification becomes 
/ ^ 0,e ^ 4/5,cri,2 2/15, '01,2 4/3. For the case of two 
chains we also list the topological symmetry sector, i.e. the eigen- 
value of the topological symmetry operator, for all energies. 

conformal weights of this state are 

hL = Af^ + + Aji + (13) 

hR = 0^/:^l + + A^:j (14) 

respectively, where Af^ = [1 + j2(j2 + l)/(^ + 1) - + 
l)/{k + 2)] is the conformal weight of the primary field 
^^^^^ su(2)^_-L X su(2)^/su(2)^. Con- 
sidering for simplicity k = odd, we can choose ji and j[ 
to run over integer values 0,l,...,(/c — l)/2 and j2, ^2 
over values 0, 1/2, 1, .., {k — l)/2. We are only interested 
in fields with Hl = Jir, so that the scaling dimension is 

We now compare the above analytical spectrum with nu- 
merical spectra for decoupled chains. Finite size scaling of 
the low-energy spectrum (similar to the procedure employed 
in previous Section) enables us to assign conformal weights 
to each energy level, analogous to the case of the 2-leg su(2)3 
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FIG. 14: (color online) Bond-energy correlations in a 4 x 8 lad- 
der with open boundary conditions (OBC) along the rungs (i.e. with 
two inner and two outer 8 -site legs) as a function of the distance be- 
tween two (parallel) bonds along the leg direction or across the lad- 
der. The bonds are oriented simultaneously along the rungs (rung- 
rung correlator) or along the legs (leg-leg correlator). The discon- 
nected part has been subtracted and the data are normalized w.r.t. the 



Jieg > i.e. in the c 

couplings, |Jrung|/Jleg 
|Jrung|/Jleg -4.38 ((9 = 



- 7/10 phase of Fig. nTTa). a) Isotropic 
= 1 (^ = 77r/4); b) strong rung coupling, 
ll7r/7). 
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FIG. 13: (color online) Low-energy excitations of 4-leg ladders with 
open boundary conditions. The ground- state energy (K=0) is used 
as energy reference. Ferromagnetic (antiferromagnetic) rung (leg) 
coupling are considered i.e. giving rise to the c = 7/10 phase of 
[1 a): moderately strong rung coupling, | Jrung | / -^leg | — 1-73 
57r/3), and strong rung coupling, | Jrung|/Jieg — 4.38 {6 = 
ll7r/7), are shown. (a,c) Low-energy spectra of a 4 x 8 ladder versus 
momentum along the ladder. (b,d) Finite size scalings of the low- 
energy levels for the two cases shown in a) and c). Fits to c=7/10 
CFT invariant spectra are provided: expected levels corresponding 
to primary (secondary) fields of the CFT are shown by red boxes 
(blue circles). The overall energy scale of the CFT is set by adjusting 
the position of the = 0, 2h = 1/5 state to the corresponding 
energy level of the 4 x 8 cluster. 



ladder in Fig. [12] As before, an accurate fit of the ground- 
state (GS) energy per site versus 1/L^ fixes the overall energy 
scale. The (allowed) combinations 2hL = 2hi + 2/i2 corre- 
sponding to the sum of two conformal weights, each arising 
from a single edge state (compare Eq. [13]), can then be read 
off from the slopes of the lowest excited states versus 1/L^. 
These numerical results as well as those for three decoupled 
legs (scaling not shown here), and for both ferro and antifer- 
romagnetic (intra-leg) couplings, are summarized in Table [II] 
The quantum numbers (and degeneracies) obtained numeri- 
cally are in perfect agreement with those obtained from the 
above analytical analysis. 



V. EFFECTS OF BOUNDARY CONDITIONS 

A characteristic feature of a topological phase isjhat it is 
sensitive to the topology of the underlying manifold^^, which 
is reflected in a non-trivial ground- state degeneracy and the 
occurrence of gapless edge modes for open boundaries. In 
this Section, we will investigate the sensitivity of the anyon 
ladder systems to these latter effects of changing boundary 
conditions. So far we considered anyon ladder systems with 
open boundary conditions along the rung direction and peri- 
odic boundary conditions along the leg directions, resulting 
in the topology of an annulus. Following the arguments in 
Refs. 15 20 we interpret the observation of gapless states in 
the energy spectrum as the appearance of gapless edge modes 
at the open boundary conditions. As a consequence, we ex- 
pect the energy spectrum to gap out as we remove the gapless 
edge states by gluing together the open boundaries of the an- 



nulus to yield a torus geometry. As detailed in Section [irD] this 
topology change is accomplished by adding a rung coupling 
between the two outer legs of the ladder (and introducing the 
correct Dehn twist). 

As an example system we consider a four-leg ladder with 
ferromagnetic rung coupling Jrung < 0. In the case of open 
boundary conditions along the rung direction, we find clear 
signatures for gapless edge modes at these open boundaries: 
First, the energy spectrum is gapless in the thermodynamic 



limit as shown in Figs. [13] The energy eigenvalues again agree 
well with the expected conformal weights, both in the strong 



and intermediate coupling regimes shown in Figs. 13:) and 
Figs. [T3k), respectively. Second, correlations of the bond- 
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FIG. 15: (color online) Low-energy excitations of 4-leg ladders with 
periodic boundary conditions. The ground state energy (K=0) is used 
as energy reference. (a,c) Spectra of a 4 x 8 ladder for values of cor- 
responding to ferromagnetic (antiferromagnetic) rung (leg) coupling. 
Isotropic couplings, |Jrung|/Jieg = 1 (^ = 77r/4), and moderately 
strong rung coupling, | Jrung|/Jieg ^ 1.73 (^ = Stt/S), are shown. 
(b,d) Finite size scalings of the low-energy levels for the two cases 
shown in a) and c) revealing 3 -fold degenerate ground states and a 
finite gap. 



energy operator decrease significantly slower between two 
bonds located on the two outer legs than on the inner legs as 



shown in Fig. 14 In addition, the bond-bond energy correla- 
tions on the rungs in both the two outer rows and the inner row 
of the 4-leg ladder decay more rapidly that their leg counter- 
parts. We interpret these differences as evidence for a gapless 
edge mode being located at the open boundaries of the ladder 
system and the presence of a gap in the bulk. 

The occurrence of the bulk gap of this anyon ladder sys- 
tem becomes even more evident when we consider the energy 
spectrum as we close the open boundary conditions, thereby 
removing the gapless (edge) modes. In Fig. [15] we show such 
clearly gapped energy spectra for periodic boundary condi- 
tions (in the intermediate coupling regime = bn/S and 
= 77r/4). This observation should be contrasted with our 
results for open boundary conditions and the same coupling 
parameters: as shown in Figs. [13^) and [13})), for the same 
coupling parameter ^ = 57r/3 the energy spectrum in the case 
of open boundary conditions nicely matches the gapless spec- 
trum of a conformal field theory. 
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FIG. 16: (color online) Low-energy excitations of 4 x 6 ladders 
with ferromagnetic rung and leg couplings. The ground state energy 
(K=0) is used as energy reference, (a) Open boundary condition (in 
the rung direction). The data are obtained for strong rung coupling, 
Jrung/ Jieg 4.38 (0 = IOtt/T). A fit to a c=4/5 CFT invariant 
spectrum is shown: expected levels corresponding to primary (sec- 
ondary) fields of the CFT are represented by red boxes (blue cir- 
cles) and the overall energy scale is set by adjusting the position of 
the K=0, 2h=4/5 energy level, (b) Periodic boundary condition (in 
the rung direction). The data are obtained for isotropic couplings, 
Jvung/ Jieg = 1 (^ = 57r/4). A large gap is seen above three quasi- 
degenerate levels. 



Furthermore, the gapped energy spectrum for periodic 



boundary conditions, as illustrated in Fig. 15 also reveals the 
occurrence of an unusual, non-trivial ground-state degener- 
acy for the anyonic ladder system. For example, in the case 
of ferromagnetic rung coupling Jrung < and antiferromag- 
netic leg coupling Jieg > 0, we observe three ground stated 
(one at momentum and two at momentum tt) separated from 
the rest of the energy spectrum by a gap of order 0(1) in the 
exchange coupling strength, which become degenerate in the 
thermodynamic limit. Evidence for the latter is provided in 



the finite-size scaling plots of Fig.[T5p) and d). It is important 
to notice that such a ground state degeneracy is not due to a 
spontaneous dimerization along the ladder direction (or to any 
other spontaneous translation symmetry breaking). Indeed, in 
the case of a spontaneous dimerization, the expected ground 
state degeneracy would be a multiple of two (depending on 
whether the system breaks translational invariance along both 
ladder directions) instead of three for the anyon ladder. 

Further evidence for a uniform anyon ground state is pro- 
vided by inspection of the correlations of the energy (rung or 



leg) bond operators shown in Fig. 17 for the same periodic 
anyon ladder. While for a dimerized system (period-2) oscil- 
lations of these correlations survive at arbitrarily large sepa- 
rations between the bonds (the amplitude is the square of the 
order parameter for infinite separation), our data show in con- 
trast a rapid vanishing of those oscillations with distance. 
Similarly, we have also checked that the spectrum of a 
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FIG. 17: (color online) Bond-energy correlations in a 4 x 8 lad- 
der with periodic boundary conditions (PBC) along the rungs as a 
function of the distance (in the leg direction) between two (paral- 
lel) bonds. The bonds are oriented simultaneously along the rungs 
(rung-rung correlator) or along the legs (leg-leg correlator). The dis- 
connected part has been subtracted and the data are normalized w.r.t. 
the zero-distance auto-correlation. Data are shown for Jrung < and 



Jieg > i.e. in the c = 7/10 phase of Fig. 11 a) and for isotropic 

1(0: 



couplings, I Jrung|/Jleg 



77r/4). 



4x6 ladder with both ferromagnetic rung and leg couplings, 
-^mng < and Jieg < 0, is fully consistent with (i) a c = 4/5 
CFT invariant spectrum when open boundaries are used, most 



evidently seen for strong rung coupling in Fig. [16] a) and (ii) 
a gapped spectrum and a three-fold degenerate ground state 
with momenta and ±27r/3 is found when using periodic 
boundary conditions (i.e. removing the edges), most evidently 



seen for isotropic couplings in Fig.[T6]b). Again, this degen- 
eracy is not connected to translation symmetry breaking but 
rather is a signature of a new (uniform) topological liquid.^^ 
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TABLE III: GS angular momentum (i.e. total spin) of a single /^rro- 
magnetic su(2)fc rung as a function of its length W (first line) and k 
(first column). The last line gives the GS total spin of ferromagnetic 
Heisenberg SU(2) open chains of same lengths. When small Jieg is 
switched on, the corresponding VK-leg ladders can be mapped onto 
single gapless effective chains (at low energies) except for the cases 
(i) marked by boxes (Haldane effective chains) or (ii) marked by "0" 
(i.e. when the GS of a single rung is a singlet). The vertical bars || 
indicate the range W < k/2 (see text). 



mod k, the total angular momentum on a rung is j = and 
we find gapped phases around this strong rung-coupling limit. 
If the ladder width W is not a multiple of the level k, i.e. 

/ mod k, then we still expect gapped phases if the 
total angular momentum on a rung is an integer (thus giving 
rise to generalized Haldane phases ^^). Similarly, we expect 
that gapless phases are found for a total angular momentum 
on a rung becoming a half-integer (and W not a multiple of 
k). These results are summarized in Table [III] 

This scenario also matches nicely the well-known behav- 
ior of ordinary SU(2) ladder models, which we recover when 
taking the limit oik ^ oo for the anyonic theories. 



VII. APPROACHING THE 2D LIMIT 



VI. SU(2)fe GENERALIZATIONS 

We now turn to the question of how the characteristic fea- 
tures of the VK-leg ladder models found for the Fibonacci the- 
ory su(2)3 are generalized when considering su(2)/c theories 
with /c > 3. All these theories allow to define ladder models 
built out of generalized angular momenta j = 1/2, similar to 
the description given in Sec. |ll] For these more general the- 
ories we can identify angular momentum j with angular mo- 
mentum k/2 — j with the highest possible allowed angular 
momentum thus becoming {k — l)/4 when k is odd. 

Following the same route as taken for the Fibonacci the- 
ory, we can access most features of their respective phase di- 
agrams by considering the strong rung-coupling limit as pre- 



sented in Sec. [Ill] In particular, such an approach reveals the 
appearance of gapped and gapless phases as a function of the 
ladder width W and the level k. For antiferromagnetic rung- 
coupling Jrung > we find that the odd/even effect of su(2)3 
occurs for all level k. On the other hand, for ferromagnetic 
rung-coupling Jrung < a more refined picture emerges: If 
the ladder width 1^ is a multiple of the level k, i.e. W = 



We conclude with a perspective on how to connect the re- 
sults obtained here for I^-leg ladders to the thermodynamic 
limit of two-dimensional lattice configurations of non-Abelian 
anyons. The strong rung-coupling limit, which was useful to 
discuss the phases of I^-leg ladders, turns out to be of little 
help in understanding this 2D limit. The reason is that the gap 
of an isolated rungs vanishes sls 1/W with increasing width, 
which restricts the applicability of the perturbative argument 
around the strong rung-couling limit to a regime of couplings 
Jieg/ Jrung < 0{1/W), wMch also vanishcs as 1^ ^ oo. 

Instead we consider the following general symmetry ar- 
gument: In contrast to their ordinary SU(2) counterpart, the 
su(2)/e anyonic theories lack a built-in continuous symmetry. 
In the assumed absence of an emergent continuous symmetry 
this reduces their ability to undergo a spontaneous symmetry 
breaking transition - such as, in two dimensions, the forma- 
tion of a Neel state and its gapless Goldstone mode for ordi- 
nary SU(2) quantum magnets. Therefore one is naturally led 
to expect gapped quantum ground states, such as topological 
quantum liquids, in these anyonic systems. This raises the 
question of how these two seemingly disjunct scenarios for 
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SU(2) and su(2)k can be reconciled when taking the /c ^ oo 
limit of the anyonic theories. Noting that the deformation of 
SU(2) used to describe the anyonic systems expHcitly breaks 
time reversal symmetry we can think of 1 //c as the strength 
of a symmetry breaking field. As such we expect the bulk 
gap of the 2D anyonic quantum ground state to close as one 
approaches the SU(2) limit, thereby smoothly connecting the 
topological quantum liquids to the Neel state. 

The formation of a gapped bulk liquid in the thermody- 
namic limit is further backed by the 'liquids picture' pre- 
sented in Ref. |20l There we have argued that the interac- 
tions between a set of non-Abelian anyons arranged on a two- 
dimensional lattice gives rise to the nucleation of a new bulk- 
gapped (i.e. topological) quantum liquid within the 'parent 
liquid' of which the anyons are excitations of. At the spa- 
tial interface between these two distinct, bulk-gapped phases 
gapless edge modes will form whose precise character can be 
identified from the gapless modes of one-dimensional chains 
of anyons which in turn allows for an identification of the 
newly formed two-dimensional bulk-gapped liquid^^. For the 
case that both the rung and leg couplings are 'antiferromag- 
netic', i.e. Jrung > and Jieg > 0, this liquid is described 
by a su(2)fe_i xsu(2)i Chern-Simons theory. On the other 



hand, if both couplings are 'ferromagnetic', i.e. Jrung < and 
Jieg < 0, then this liquid is described a [/(I) Chem-Simons 
theory. The case of mixed coupling signs remains open. 

Anyonic generalizations of quantum magnets in the spirit 
of the work presented here can discussed in analogous fash- 
ion for other anyonic theories (tensor categories) and for other 
two-dimensional lattice geometries and interactions. We ex- 
pect this to be a fruitful and broad field of research at the inter- 
face of quantum magnetism and topological states of matter. 
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